Viscosity of liquid "^He and quantum of circulation: Why and how are they related? 
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The relationship between the apparently unrelated physical quantities - kinematic viscosity of 
liquid *He, v, and quantum of circulation, k = 27r/i/m4, where h is the Planck constant and m4 
denotes the mass of the '*He atom - is examined in the vicinity of the superfluid transition occurring 
due to Bose-Einstein condensation. A model is developed, leading to the surprisingly simple relation 
v ~ k/6. We critically examine the available experimental data for ^He relevant to this simple 
relation and predict the kinematic viscosity for the stretched liquid *He along the A-line at negative 
pressures. 
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In this Letter, following Lars Onsager [1] who was the 
first to raise this issue, we examine a physical relation 
between two quantities of apparently completely differ- 
ent physical meaning: the kinematic viscosity of normal 
viscous fluids, responsible for their classical mechan- 
ical property - friction between two fluid layers moving 
with different velocities; and the quantum of circulation. 
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where Ms is the mass of the superfluid particle (such 
as the mass of ''He atom or two '^He atoms) describing 
quantization of velocity circulation in integer multiples 
of K. Both quantities v and k have the same dimension, 
[v] = [k] — cm^s~^. This rises the immediate question 
whether or not v and k are related in fluids experienc- 
ing transition to superfluidity such as ^He (at A-point 
T = Tx), ■^He or Bose-Einstein condensates (BEC) of 
cold clouds of alkali atoms, and if so, why and how? It 
should be noted that is generally temperature and pres- 
sure dependent, while for a given medium k is a physical 
constant. We therefore bound the temperature region 
close to the physically relevant temperature T* (which in 
the particular case of "'He is = T\). Moreover, in order 
to guarantee the conventional viscous fluid behavior, we 
restrict the temperature region hy T > Tx and, in order 
to neglect compressibility effects, we do not consider high 
external pressures far above the saturated vapor curve. 

For classical hydrodynamic community the question of 
a possible relationship between v and k seems irrelevant. 
Classical hydrodynamics does not involve the Planck con- 
stant, h, and quantization of any relevant quantity does 
not take place. At the same time, fluid viscosity, v, does 
not concern most researches in quantum physics. There- 
fore we should address this question to "general" theo- 
reticians and to the "quantum" superfluid community. In 
particular, this question is of high interest for researchers 
in the field of quantum turbulence [2-4] , where the prob- 
lem of turbulent energy dissipation in superfluids is a 
hot topic [6-16]. It is usually discussed in terms of effec- 
tive kinematic viscosity in units of jk, where the pref- 



actor 7 is expressed via phenomenologically introduced 
parameters such as mutual friction coefficients or the Kol- 
mogorov constant. At the same time, as far as we know, 
the relation between v and k has not yet been seriously 
addressed. 

At first glance, being of very different physical origin, 
v and K can be completely different. It is known (see 
Fig. 1) from experiments that in the normal liquid He-I 
first increases slightly along the saturation vapor pressure 
(SVP) curve from about 2.59 x 10^''' cm^s^'^ at its normal 
boiling point close to 4.2 K and then drops down to about 
1.67 X 10"'' cm^s~^ just above Tx ~ 2.17 K, i.e., v varies 
by about 30% and is of order 0.1 k [19, 20]. The super- 
fiuid community tends to consider this experimental fact 
as accidental [2] , often having in mind the very different 
situation in fcrmionic ■^He: Along the SVP the kinematic 
viscosity decreases slightly, but its lowest value is still 
about an order of magnitude bigger than that of He I in 
this range (see Fig. 1), but with decreasing temperature 
v increases and above the superfiuid transition in the mK 
range it behaves as a thick Fermi liquid (for which the 
Landau Fermi liquid theory [21] predicts v oc T^^) with 
the (temperature and pressure dependent) kinematic vis- 
cosity comparable with that of air or machine oil, thus 
i^/k ^ I close to the critical temperature, Tc. 

"General" theoretical physicists think differently: in 
the discussed temperature region both v and k (and Tx 
as well) can be linked uniquely using atomic units, such 
as the Planck constant h, the atomic mass M or the Bohr 
radius. For v and n this gives [in agreement with Eq. (1)] 
V ^ K ^ h/M, or ly/K ~ 1, i.e. "of the order of unity, at 
least in Bose-fluids (such as ''He), while in Fermi- fiuids 
(such as ^He) this is not necessarily the case [22]". 

This answer, for *He in particular, cannot completely 
satisfy "practically oriented researchers" for whom num- 
bers such as 47r w 12 or l/27r « 1/6 (which can often 
be found in physical equations) cannot be considered as 
having the same order of magnitude, although from the 
general theoretical viewpoint this is considered to be the 
case. 
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FIG. 1. (Color online) The temperature dependence of the 
kinematic viscosity of liquid *He (based on total fluid den- 
sity) and liquid ^He (calculated using the dynamic viscosity 
as measured and fitted by Black, Hall and Thompson [17] and 
density values as analyzed and fitted by Huang, Chen, Li and 
Arp [18]) along their SVP curves. Note that the minimum v 
value in normal liquid ''He exceeds that in normal liquid *He 
by a order of magnitude. 

In this Letter, we use very simple physical arguments 
instead of just dimensional reasoning in order to find 
some more quantitative relationship between v and k 
rather then just z/ ~ k, especially for the particular case 
of the most common quantum fluid - liquid ^He. 

The model. We start with the Boltzmann kinetic 
theory for weakly non-ideal gases, in frame of which v 
can be estimated as: 

V K, -v \ . (2a) 

Here A denotes the mean-free path and v is the thermal 
velocity, at temperature T >T^. given by 

Mv^ 3 , , 

^=2fcBT, (2b) 

where is the Boltzmann constant and M denotes the 
atomic mass. Equations (2a) and (2b) give 

Kr).y^fA, (2c) 

Let us consider liquid *He and substitute for T its 
particular value T» = Ta, estimated by the celebrated 
London-1938 formula for T^, the temperature of Bose 
condensation of ideal gas [23]: 

/ n 27rfi2 ;i2„2/3 



where C,{x) is the Riemann zeta function; C(3/2) ~ 2.6124 
and n denotes the number of bosons in a unit volume, 
n = N/V. Based on the experimental data available 
at that time, London evaluated from his Eq. (3) T^^ ~ 
3.1 K. As it was pointed out many times, taking into 
account the simplicity of the ideal Bose gas model, this 
value is strikingly close to the experimentally observed 
Tx « 2.2 K. Having in mind that the temperature enters 
Eq. (2c) under the square root, one estimates the possible 
error of this replacement as — 1 ~ 0.2, i.e., within 

about 20%. Optimists may say that this disagreement 
(up to 20%), originating from interatomic interactions, 
may also be considered as an inaccuracy estimate for our 
simple model of the kinematic viscosity. 

The next step is to estimate the mean-free path. A, 
which in weakly interacting gas exceeds the mean inter- 
atomic distance ^ = n~^/^ considerably. Due to physical 
reasons (that lie outside the scope of this discussion), 
A « £ in the densely packed liquid '^He at temperatures 
T <T\, with a good accuracy (much better than 20%). 
In this way, we implicitly account for strong interaction 
between liquid-helium atoms, thus exploiting the appli- 
cability limit of the Boltzmann Eq. (2a) . 

Now, Eqs. (2c), (3) and (1) allow one to estimate the 
kinematic viscosity of ^He just above its superfluid tran- 
sition as 

,^ - /3.31 h /3.31 K K 

i^siTx) ~ \ = A/ w0.167kw-. (4) 

•^'^ V 3 Af V 3 27r 6 ^ 

Let us compare this simple model prediction with the 
available experimental data. In fact, the thermodynam- 
ical and transport properties of cryogenic ^He are very 
well known. Along the saturated vapor pressure (SVP) 
curve (see Fig. 1) they have been tabulated by Donnelly 
and Barenghi [19]; we use their value of Tx in Table I. 
Equilibrium thermodynamical properties of ^He can also 
be accurately calculated, thanks to Arp and McCarty, 
using the program HEPAK [20]; which we have used to 
evaluate values of kinematic viscosity, z/, at all positive 
pressures given in Table I. 

Having in mind the approximative character of our rea- 
soning, the accuracy of which we optimistically estimated 
as 20%, we have to consider the surprisingly exact agree- 
ment between the estimate (4) and the experimental data 
at saturated vapor pressure (SVP) given in Table I as un- 
expectedly good and accidental. In any case, even 20% 
of inaccuracy is a much better result than the simple di- 
mensional estimate v{Tx) ~ k. 

What can we say in brief about London's estimate 
of Txj used in our model? As we mentioned, he found 
Tl about 3.1 K in the approximation of ideal Bose gas 
which corresponds to the low density limit, when in- 
teratomic interactions can be neglected. On the other 
hand, it follows from the influential theoretical work of 
Bogolyubov on weakly non-ideal Bose gas [24] that inter- 
actions are essential for the phenomenon of superfluidity. 
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2.095 


0.1580 
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1.646 


2.157 


0.1492 


1.72 X 10" 


4 


0.173 


SVP 


2.172 


0.1461 


1.67 X 10" 
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0.168 


spinodal 


~ 2.2 


~ 0.116 


~ 1.6 X 10 


-4 


0.167 



0.35 



0.30 



TABLE I. Selected properties of liquid ^He along the A-line. 
For details, see the text. 



as for ideal Bose gas with quadratic dispersion relation 
the well-known Landau criterion [25] gives zero critical 
velocity! From a practical point of view, the lowest pos- 
sible density corresponds to metastable liquid along the 
nucleation line close to the negative spinodal pressure, 
at which T\ « 2.2 K, which is again quite close to the 
London estimate. The approximation of ideal gas grad- 
ually fails with increasing density. It is therefore natural 
to think that the experimentally observed T\ will de- 
viate from the predicted with increasing density of 
liquid '*He, i.e., with externally applied pressure, as the 
real system of liquid "^He deviates more and more from 
the ideal model of non-interacting Bose gas. This is in- 
deed the case: the experimentally observed values are not 
proportional to the density in power 2/3 as predicted by 
London, moreover, with increasing pressure they even 
decrease, making the slope dp/dT\{p) negative. To illus- 
trate this, we give the observed values of T\ at several el- 
evated pressures, p, based on the accurate measurements 
of Ahlers [26] in Table L 

Fig. 2 plots the dimensionless ratio v/ n for liquid ''He 
at temperatures where the superfluid transition is ob- 
served at various applied pressures. For the circulation 
quantum we use k = 0.997 x 10"'^ cm^s~^. We see that 
y / K, behaves in a regular way and with decreasing the 
external pressure approaches the value = 0.167 given by 
Eq. (4) . It appears that our simple model works surpris- 
ingly well at low external pressures. This fact enables us 
to predict the value of kinematic viscosity for metastable 
liquid ^He that could exist at lower pressure than the 
equilibrium SVP and even at negative pressures down to 
the nucleation line, where cavitation occurs (for review 
on cavitation in liquid helium, see Balibar [27] and ref- 
erences therein). The location of the A-line at negative 
pressures has been calculated by Maris and Edwards [28] . 
The maximum value of Tx = 2.205 K is reached at a neg- 
ative pressure of about -5.7 bar, and the A-line meets 
the spinodal line at a negative pressure of -7.3 bar and 
a temperature of about 2.18 K. Along this A-line in the 
metastable liquid, the kinematic viscosity value ought to 
be close to 0.167 n. 
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FIG. 2. (Color online) The ratio v/k plotted versus exter- 
nal pressure using values of Table I compared with the value 
predicted by Eq. (4) shown as dotted (blue) line. Negative 
pressures down to spinodal pressure, where metastable liq- 
uid could exist (red, left of the vertical arrow marking SVP), 
are included. Our approach allows to estimate the kinematic 
viscosity (see the red thick solid line) for the metastable liq- 
uid along the A-line, which ought to extend down to the nu- 
cleation line. Its position in the pressure-temperature phase 
diagram was calculated [28], but not as yet experimentally 
accurately determined. 



Bosonic substances. Our approach would allow, in 
principle, the prediction of kinematic viscosity for other 
liquid bosonic substances in which the phenomenon of 
Bose-condensation could occur. However, clouds of spin 
polarized hydrogen or alkali atoms at conditions when 
BEG happens represent very dilute gases and our above 
estimate of mean free path cannot be used here. For 
fermionic substances the situation is more compli- 
cated [29]. 

We are fully aware that our approach could be right- 
fully criticized, both from the standpoint of general theo- 
retical physics as well as from the point of view of classical 
fluid dynamics. Nevertheless, having in mind that the- 
ory of superfluid turbulence (especially problems of tur- 
bulent energy and/or vortex line density decay involving 
phenomenologically introduced quantity - effective kine- 
matic viscosity) is still in its infancy, we hope that our 
approach can be considered as its useful building block. 

To conclude, we have considered the physical relation 
between the two apparently independent physical quan- 
tities - kinematic viscosity, i^, versus quantum of circula- 
tion, K - and developed a simple model predicting that for 
bosonic liquids at the temperature of Bose condensation 
they are related as u ~ k/6. We have critically examined 
the experimental data on the transition to superfluid- 
ity at various external pressures for the most common 
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bosonic liquid - ^He - and found that in the relevant 
region of parameters (T ~ Tx and at low applied pres- 
sures) the agreement is » 10^ times better then one may 
expect from dimensional analysis, which is usually accu- 
rate within about an order of magnitude. We also pre- 
dict that this simple relation ought to hold even better 
for the mctastablc stretched liquid along the A line at 
negative pressures. One can estimate the probability of 
this agreement, ~ 10~^, as totally accidental. This raises 
many questions, for example: (i) Why does it happen? 
(ii) Why do ideal-gas considerations work so well for the 
real liquid? (iii) Could it be thanks to cancelation of 
effects of some as yet not known factors, which are defi- 
nitely ignored in our suggested approach? (iv) Are these 
cancelations accidental or are there some deep physical 
reasons yet to be discovered? In other words, our Let- 
ter asks more general questions than it provides clear 
answers. This, we hope, will trigger new discoveries in 
fields of physics that are considered old and apparently 
well established. 
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[29] In case of fcrmionic substances wo need to consider two 
characteristic teinporaturcs of purely quantum mechan- 
ical origin: the Fermi temperature Tp and the criti- 
cal temperature Tc at which the superfluid transition 
occurs (usually Tc <C Tp). The crossover from classi- 
cal behavior at high temperatures to quantum mechan- 
ical behavior at low temperatures is demarcated by Tp, 
as the subsequent superfluid (superconducting) transi- 
tion (in conventional BCS sense) is a purely quantum 
mechanical effect. Let us consider ^He and substitute 
for T its particular value T, = Tp, estimated by the 
well known formula for the Fermi energy of the ideal 
Fermi gas: = /cbTp = h^ {Q-w'^nY'^ /{2M) . In anal- 
ogy with the approach outlined above we formally ar- 
rive at z^(Tf) fs k/ (67r*)^^'' fs 0.346 /t , where we took 
into account that the mass of the superfluid particle is 
Ms = 2M; for ^Ho k 0.66 x 10^^ cm^s"\ By com- 
parison with a minimum in Fig. 1 we see that for the 
fermionic "'He the predicted value of kinematic viscosity 
is about an order of magnitude off. This disagreement is 
perhaps not too surprising, however, because superflu- 
idity occurs at the transition temperature three orders 
of magnitude lower, where quantum mechanical effects 
arc dominant and our starting assumption that the kine- 
matic viscosity can be estimated from the Boltzmann ki- 
netic theory for ideal gases is not justified. Instead, as 
it has been shown by Landau in frame of his Fermi liq- 
uid theory (see, e.g., [21]), in the temperature region of 
interest kinematic viscosity behaves as oc and as 
the temperature approaches Tc our simple model is not 
applicable. 



